Absence of localized-spin magnetism in the narrow-gap semiconductor FeSb 2 . 
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We report the inelastic neutron scattering measurements aimed at investigating the origin of 
temperature- induced paramagnetism in narrow-gap semiconductor FeSt>2. We find that inelastic 
response for energies up to 60 meV and at temperatures « 4.2 K, w 300 K and ~ 550 K is 
essentially consistent with the scattering by lattice phonon excitations. We observe no evidence for 
a well-defined magnetic peak corresponding to the excitation from the non-magnetic S = singlet 
ground state to a state of magnetic multiplet in the localized spin picture. Our data establish the 
quantitative limit of Sf// < 0.25 on the fluctuating local spin. However, a broad magnetic scattering 
continuum in the 15 meV to 35 meV energy range is not ruled out by our data. Our findings make 
description in terms of the localized Fe spins unlikely and suggest that paramagnetic susceptibility 
of itinerant electrons is at the origin of the temperature- induced magnetism in FeSb2. 

PACS numbers: 71.28+d, 75.40.Cx, 75.40.Gb, 75.50. Ee 
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I. INTRODUCTION 

A tight balance between strong covalent hybridiza- 
tion, the tendency to band delocalization, strong elec- 
tronic correlation and the crystal field potential leads to 
a rich variety of unusual electronic states in the inter- 
metallic pnictides and chalcogenides TM X , where T is 
the transition metal (Mn, Fe, Co, Ni, ...), M is a met- 
alloid (Sb, As, Te, ...), and 0.5 < x < 3^Among 
them are semi-metallic and semiconducting compounds 
with rich variety of magnetic properties and magnetic 
structures^— half-metallic ferromagnetS ) 10 ' 1 1 materials 
showing giant magnetocaloric effect ^2. unusually high 
thermopowerfi 3 - - — colossal magneto-resistance and giant 
carrier mobility Many of the unique electronic prop- 
erties which are of great current and future technological 
interest are intertwined with the emerging magnetism in 
these materials. Recent discovery of the new variety of 
the high-temperature superconductivity in layered tran- 
sition metal pnictides, which is also closely related with 
magnetism,— ~— has generated new surge of interest in 
the physical properties of these compounds. 

In defiance of the naive expectation that alloys com- 
posed of metallic constituents should be metallic, many 
pnictides and chalcogenides are actually semiconducting, 
or even insulating. Non-metallic and half-metallic behav- 
iors result primarily from the hybridization gaps imposed 
by strong covalent bonding between the transition metal 
3d and the metalloid p (or, more precisely, sp) orbitals, 
which dominate the electronic band structure near the 
Fermi level. These are also further enhanced by strong 
correlation of the narrow-band d- elect ronsi 2 ^ 2 ^— 

Iron diantimonide FeSb2 is a narrow gap semiconduc- 
tor representative of the TM2 pnictide family with a 
number of interesting and unusual behaviors* ^ 18 ' 29 " — 
Perhaps most outstanding are the temperature-induced 
paramagnetism similar to that found in iron silicide, 
FeSi^ 3 - and strongly anisotropic electrical transport 



properties Below room temperature and down to 
about 40 K the conductivity of single crystal FeSb2 shows 
metallic behavior along one of the crystallographic direc- 
tions (the resistivity for current along the c-axis decreases 
with the decreasing temperature) and semiconducting ac- 
tivated behaviors along two other, a and b axes 3 - 2 -. Such 
behavior suggests anisotropic band structure, which may 
results from a particular network of hybridized orbitals, 
a flavor of the orbital order. Then, either there is a one- 
dimensional metal-insulator transition changing the o 
axis transport around 40 K, or a sequence of semiconduc- 
tor cross-overs corresponding to anisotropic band gaps, 
which vary from about 4 meV for transport along the c 
axis to about 30 meV perpendicular to it. 

As conductivity of FeSb2 increases with the increas- 
ing temperature, an excess Shottky-like electronic heat 
capacity 3 ^ and an exponentially-activated paramagnetic 
susceptibility 2 ^ 31 ^ 35 also develop. The origin of these 
and other intriguing properties, such as the giant ther- 
mopower, colossal magneto-resistance and high carrier 
mobilities, lies in the physical nature of the lowest 
unoccupied and the highest occupied electronic states 
near the Fermi level. Currently, the structure and the 
composition of these states are unclear. In particu- 
lar, the temperature-induced paramagnetic susceptibil- 
ity could be analyzed either in the framework of the 
ionic, localized-spin model, or as a Pauli paramagnetism 
of itinerant electrons belonging to two bands separated 
by a narrow gap containing the Fermi energy, Ep. 31 
In the former, localized-spin picture, it involves ther- 
mally activated spin state transition between the non- 
magnetic S = and magnetic S=l states of the crystal- 
field split multiplet of Fe 3c? electronic levels, correspond- 
ing to semi-localized narrow bands^ 9 - In this scenario, the 
corresponding crystal-field spin excitation should be ob- 
servable in the inelastic magnetic neutron scattering ex- 
periment. Search for this excitation aimed at clarifying 
the nature (localized vs itinerant) of the temperature- 
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induced magnetism and thus the lowest unoccupied elec- 
tronic states in FeSb 2 , was the goal of the present study. 

Fitting the temperature-dependent magnetic suscepti- 
bility x{T) to a thermally induced spin-state transition in 
the localized spin model yields a spin gap A x w 47 meV 
and an effective magnetic moment of about \.2\xb (Bohr 
magnetons) for the magnetic high-spin stated The cor- 
responding singlet-triplet spin excitation would therefore 
be expected to produce a peak of significant intensity in 
the FeSb2 magnetic neutron scattering cross-section at 
an energy « A x . The intensity of such an excitation is 
expected to show characteristic temperature dependence. 
Namely, it should decrease with the increasing temper- 
ature, as the transition probability gets depleted with 
the increasing thermal population of excited magnetic 
states.— We have conducted an extensive search for such 
magnetic spin excitation for energies up to 60 meV and 
at temperatures 10 K, 300 K and 550 K using inelastic 
thermal neutron scattering, which is reported here. 

The paper is organized as follows. Section 2 contains 
description of the FeSb2 crystal structure and discussion 
of its relation with the electronic properties of this mate- 
rial (the detailed qualitative discussion of band hybridiza- 
tion aspects illustrated in Fig. Q] follows Goodenough's 
paper [2|; it is rather specialized and could be skipped 
on the first reading). Experimental details are described 
in Section 3, the resulting data and data analysis are 
presented in Section 4, which is followed by the brief 
summary of our conclusions in Section 5. Appendix A 
presents details of the quantitative comparison of mag- 
netic and phonon neutron scattering intensities. 



II. FeSb 2 CRYSTAL AND ELECTRONIC 
STRUCTURE 

Iron diantimonide FeSb2 crystalizes in a FeS2 marc- 
asite structure shown in Figure Q] (a)^ ' 29 i 31 It features 
chains of edge-sharing FeSbg octahedra running along the 
c axis, which are stacked in a somewhat corrugated body- 
centered geometry, sharing corners in a and b directions. 
The lattice is orthorhombic (space group Pnnm), with 
two formula units per unit cell and room temperature 
lattice parameters a 5.83 A, b s» 6.53 A and c « 3.2 A. 

The electronic structure of Fe 3d levels in FeSb 2 is 
determined by the hybridization with the Sb bp and 5s 
orbitals (or sp 3 hybrid orbitals, but we retain p x ,p y ,p z 
tags for simplicity) and by the crystal field of the Sb octa- 
hedral environment. At low temperature, Sb octahedra 
are slightly squashed, with two shorter ("apical", « 2.56 
A below w 100 K) and four longer ("planar", ss 2.59 A 
below ps 100 K) Fe-Sb bonds. 31 Such distortion is oppo- 
site to that found in the family of high-temperature su- 
perconducting layered perovskite cupric oxides and their 
nickel- or cobalt-based relatives, where oxygen octahedra 
are elongated and tetragonal crystal field favors d yz . zx 
and G?3 2 2_ r 2 orbitals. Schematics of the resulting band 
structure following from qualitative arguments of Ref. 2 




FIG. 1: (a) Crystal structure of FeSt>2. Atoms beyond the 
unit cell are shown to illustrate the main structural motif - 
Fe-Sb octahedra sharing edges along c direction and corners 
along a and 6. (b), (c) Schematics of the electronic band 
structure of FeSb2 resulting from the crystal field splitting 
and strong covalent hybridization, (c) illustrates the spin- 
state transition and the corresponding crystal field excitation 
in the quasi-localized narrow-band scenario. 

is illustrated in Fig. Q](b), (c). 

In discussing the local Fe 3d electronic structure, it is 
convenient to associate the long and the short bond di- 
rections with (X, Y) and Z coordinate axes, respectively. 
In this notation, the tetragonal crystal field of squashed 
octahedra splits the e g and ti g multiplets so that d x 2_ y 2 
and d xy orbitals have lower energy compared to d 3z 2_ r 2 
and d yz ^ zx orbitals, respectively. Such level hierarchy is 
opposite to that used in Ref. for constructing schematic 
band structure of marcasite iron pnictides starting from 
the ionic Fe 4+ 3d 4 configuration. This distinction, how- 
ever, is unimportant, since distortion of the octahedra 
is quite small and further decreases upon heating, es- 
sentially vanishing at T w 500 K.— Hence, the electronic 
band structure is mainly determined by covalent shifts re- 
sulting from 3d— bp hybridization and formation of bond- 
ing and antibonding states. This pushes c^z^-derived 
states down, perhaps below the non-bonding or weakly- 
bonding d xy level, so that splitting of the ti g multiplet is 
the same as for the crystal field of elongated octahedra. 
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An interesting problem arises with assigning an ionic 
state, such as 3d 4 , to Fe. With such an assumption, 
the Fermi level lies within the ti g multiplet and the 
lowest unoccupied states arise from the essentially non- 
bonding d xy orbital forming a narrow, semi-localized 
conduction band^ Such scenario could in fact explain 
the anisotropic conductance experimentally observed in 
FeSb2. 2 ^ Indeed, a band derived from weakly overlapping 
nearly non-bonding d xy orbitals could give rise to nearly 
one-dimensional conductance along the c-axis direction, 
Fig. [U Such assignment, however, disagrees with the 
recent LDA + U band structure calculations 1 27 i 28 which 
find that the Fermi level lies above t^g states, and the 
conduction band is mainly of d% z 2_ r 2 origin, as shown in 
Fig. [H(b), (c). Ionic picture also disagrees with the gen- 
eral expectation of strong 3d — bp covalent Fe-Sb bonding 
in a metal-metalloid compound. It is also possible that 
covalent shifts are even stronger than those shown in Fig. 
[T] (b), (c), so that d 3z 2_ r 2 and d x 2_ y 2 levels appear be- 
low the d xy level. While the precise level hierarchy is 
presently unclear, it does not impact the arguments dis- 
tinguishing between the itinerant and localized-spin mag- 
netism presented below. 

In the hybridization-dominated covalent band scheme 
of FeSb2 illustrated in Fig. [T](b), (c), the d a .2_ y 2-derived 
valence and d 3z 2_ r 2 -derived conduction bands agree with 
the LDA + U results of Refs. H3, H The difference 
between panels (b) and (c) is in the bandwidth of the 
valence and the conduction bands near the Fermi energy, 
Ep. In the case of Fig. [T] (b), thermally excited elec- 
trons and holes are strongly itinerant and exhibit a weak 
Pauli paramagnetism, similar to the picture proposed for 
FeSii^L^I There is no local spins and local magnetic mo- 
ments in this scenario, and only a weak and broad signal 
is expected in a magnetic neutron scattering experiment. 
This is consistent with the weak temperature-induced 
quasi-elastic paramagnetic scattering observed by neu- 
trons in FeSij2£ In the narrow-band scenario of panel (c) 
on the other hand, the strong on-site correlation domi- 
nates the electronic properties, electrons near the Fermi 
level are essentially localized, and temperature-induced 
magnetic state corresponds to a local S=l triplet. In the 
band picture language this can be viewed as a localized 
S = 1 electron-hole triplet bound state. A well-defined 
peak in magnetic neutron scattering at an energy corre- 
sponding to the local singlet-triplet spin excitation, and 
with the characteristic decrease in intensity with the in- 
creasing temperature, could then be expected.— 



III. EXPERIMENTAL PROCEDURE 

Our FcSb2 crystals were similar to those used in the 
previous studies reported in Refs. [Tt], [H, HH, and were 
grown from the excess antimony flux as described in Ref. 
|29 . Our sample was an array of four large single crystals 
with total mass m = 7.46 g mounted on an aluminum al- 
loy (A16061) sampleholder shown in Fig. [5] (a). Crystals 



were co- aligned to within 0.7° in the horizontal scatter- 
ing plane, with a lattice direction being vertical. The 
sample assembly shown in Fig. [2] (a) was mounted in the 
closed-cycle refrigerator capable of maintaining temper- 
ature in the range from about 4 K to about 600 K. In all 
measurements the (0, k, I) reciprocal lattice plane of the 
sample was kept in the horizontal scattering plane. 
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FIG. 2: (a) All-aluminum sample holder with 4-crystal FeSb2 
assembly used in our measurements. (b) Scaled to A -1 
schematics of the (Okl) reciprocal lattice zone. Grey bars 
show scan directions for maps shown in Figures [3] and [S] 



Neutron scattering measurements were performed us- 
ing the HB1 thermal neutron triple axis spectrometer 
at the HFIR Center for Neutron Scattering at the Oak 
Ridge National Laboratory (ORNL). Monochromatic in- 
cident neutrons were obtained using the (002) reflec- 
tion from vertically focussing pyrolytic graphite (PG) 
crystals, scattered neutrons were analyzed using similar 
PG(002) analyzer crystals. Two 1" thick PG filters were 
placed after the sample to suppress scattered beam con- 
tamination by higher order reflections in PG monochro- 
mator and analyzer. Neutron beam collimations were 
« 48' - 40' - 60' - 120', from reactor to detector. Two 
final scattered neutron energies were used, Ef = 14.7 
meV for the high-resolution mode, and Ef — 30.5 meV 
for the low-resolution and high-intensity mode. In the 
latter case, the total volume of the sampled phase space 
increases by a factor between 4 and 5 in our energy and 
momentum transfer range. This results in a proportion- 
ally higher sensitivity to weakly dispersive features in 
the scattering cross-section. Broad surveys of scattering 
intensity for energy transfers up to 60 meV and at tem- 
peratures ~ 4.2 K, R3 300 K and ps 550 K were performed 
using high- intensity configuration with Ef = 30.5 meV. 
Lines in Fig. [5] (b) show the corresponding constant- 
energy scans in the b* — c* reciprocal lattice plane. 
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TABLE I: Intensity, position and the full width at half maximum (FWHM) of the resolution-corrected Lorentzian fits to four 
distinct phonon groups shown in Fig. [3] FWHM < 0.1 meV is less than 10% of the instrument energy resolution and implies 
no experimentally observable width. 

Wave vector, Integral Intensity (cts-meV) Position (meV) FWHM (meV) 

Q(r. 1. u.) " Ji h h h -Ei E 2 E 3 E 4 Wi W 2 W 3 W 4 

(0,0,2.5) 133(10) 130(10) 41(6) 134(12) 16.5(2) 19.7(2) 23.8(3) 30.2(2) 1.5(4) 1.2(4) <0.1 < 0.1 

(0, 4, 0) 24(5) 14(5) 53(7) 63(7) 19.4(3) 24.0(5) 29.7(3) 34.7(3) < 0.1 < 0.1 < 0.1 < 0.1 



IV. RESULTS AND DISCUSSION 

We first searched for the well-defined crystal-field 
singlet-triplet excitation, which would be expected in the 
narrow-band localized-spin picture, by carrying out en- 
ergy scans at different wave vector transfers in the high- 
resolution configuration with Ef = 14.7 meV. Several 
such scans for wave vectors Q near b* and c* directions 
are shown in Figure [3J We observed no features which 
could be identified with the expected magnetic excitation 
in the 35 meV to 60 meV range. A slight increase in the 
background (BG) towards higher energies, which is more 
pronounced for smaller Q, is an instrumental effect asso- 
ciated with the detector vessel approaching the incident 
neutron beam. This was confirmed by background mea- 
surements with sample removed from the beam, which 
are shown by solid black circles in Figure [3J (a, b) . 

A number of peaks seen in the energy range from 
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FIG. 3: Constant-Q scans near Q = (0, 0, 2), ~ along the 
c-axis direction (a) and near Q = (0, 4, 0), « along the b- 
axis (b), measured with E f = 14.7 meV at T = 4.1(1) K. 
Black circles show background measured without the sample. 
Lines show fits to a number of resolution-corrected Lorentzian 
peaks describing phonon scattering. Horizontal bars show 
the calculated full width at half maximum (FWHM) of the 
instrument resolution,^ projected on the energy axis. 



10 to 35 meV are attributable to optic phonon modes. 
A plethora of optic phonons were observed at these 
energies by Raman and far-infrared (FIR) optical 
spectroscopy^ - — Six phonon modes at 18.7 meV, 19.0 
meV (A g symmetry), and 11.2 meV (B 2g ), 18.8 meV 
(B^g), 19.1 meV and 21.6 meV (Bi g ) were identified 
in Raman experiments , 39 i 40 Additional peaks, including 
high-energy modes, were observed in the FIR reflectivity 
measurements at 13.2 meV, 28.6 meV, 31.9 meV, 33.6 
meV (electric field of light || b) and 15.0 meV, 26.8 meV, 
32.4 meV (electric field of light _L b). Solid lines in Figure 
[3J (a, b) show resolution-corrected^ fits of our data to a 
number of Lorentzian peaks, whose parameters are listed 
in Table HI Within the instrumental resolution shown 
by horizontal bars in Figure [3J peak positions observed 
in our neutron scattering measurements agree well with 
Raman/FIR data. Significant intrinsic width of E\ and 
E2 peaks at Q = (0, 0, 2.5) reflects the fact that there are 
several distinct phonon modes within each peak. 

The phonon origin of peaks in Figure [3J is further sup- 
ported by the temperature dependence of their inten- 
sity, which increases upon heating. Such behavior is 
indeed typical of phonons, whose contribution to scat- 
tering cross-section at a wave vector Q is given by the 
dynamical correlation function S P (Q 1 E) of an oscillator, 
ft is related to the imaginary part of the oscillator dy- 
namical susceptibility Xp(Q> E) through the fluctuation- 
dissipation theorem ) 42 ' 43 

7rS p (Q,E)= X ;(Q,E)/(l-e- E / k * T ) . (1) 

In the temperature range where Xp(Q, E) of the phonon 
oscillator does not change much (damping, anharmonic- 
ities, etc, are small), the measured S P (Q,E) increases. 
This increase is due to the decrease of the detailed 
balance factor (often called Bose factor, since (1 — 
e -E/k B Tyi = i + n B {T)) in the denominator of Eq. [TJ 
and is most pronounced for E < T '. The imaginary part 
of dynamical susceptibility of a damped oscillator can 
be represented as a difference of two Lorentzian peaks 
centered at ±£b [Eo/H is the undamped oscillator fre- 
quency; this naturally satisfies the causality requirement, 
Xp(Q, E) = — Xp{Qi —E)). In cases like ours, where peak 
energy is much larger than damping, the contribution of 
the negative-energy peak can be neglected at E > T, 
which justifies using simple Lorentzian line shapes. 

Having established the phonon nature of dominant 
scattering features observed in FeSb2 at low T in scans 
with Ef = 14.7 meV, we decided to perform broad sur- 
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FIG. 4: Contour map of neutron inelastic scattering intensity 
for Q = (0, 0, 0, (a-c), Q = (0, 3.1, £), (d-f), and Q = (0, 1 + 
£,£), (g-i), at T = 5 K, 300 K and 550 K (from left to right). 
Intensity is shown in counts per monitor count corresponding 
to counting time of ~ 1.5 minutes at 4 meV and ~ 4 minutes 
at 60 meV. Ellipses show the calculated full width at half 
maximum (FWHM) instrument resolution.— 



veys of scattering using the high-intensity mode with 
E f = 30.5 meV. By studying its evolution with temper- 
ature we could then attempt to single out some evidence 
of the non-phonon magnetic scattering. To this end, we 
have measured scattering along (0,0, £), (0,3.1,£), and 
(0, 1 + directions shown in Fig. [2] (b), for energies 
up to 40 meV, or 60 meV. Corresponding color contour 
maps of the measured inelastic scattering intensity are 
shown in the three rows, (a-c), (d-f) and (g-i), of Figure 
|4j The data was collected in a set of constant-energy 
scans with step dQ ss 0.075 A, taken every 2 meV. Three 
columns in the figure show intensities measured at three 
different temperatures, T = 6 K, 300 K and 550 K. El- 
lipses in panels (a), (d) and (g) illustrate the calculated 
full width at half maximum (FWHM) of the instrument 
resolution function projected along the energy axis4i 

Similar to the Ef — 14.7 meV data of Fig. [3J the scat- 
tering in Fig. 2] is also dominated by phonons and in- 
creases strongly with the increasing temperature. There 
is no evidence for any scattering above the BG level in 
the 35 meV to 60 meV range. Both optic phonons and 
an acoustic phonon emerging from the (0,0,2) Bragg 
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FIG. 5: Imaginary part of the dynamical susceptibility ob- 
tained by applying the detailed balance factor of Eq. [T]to the 
data shown in Figure [4] 

peak are seen in (0,0, £) maps of Fig. [4] (a-c), while 
for (0, 3.1, £) and (0, 1 + £, £) the scattering is dominated 
by optic phonon modes. The imaginary part of the dy- 
namical susceptibility obtained from the data of Fig. @] 
upon subtracting the measured (Q, E)- independent BG 
of 4 counts/min and using Eq. Q] is shown in the cor- 
responding panels of Fig. [5j It quantifies the system's 
oscillator response, free of the extra temperature depen- 
dence of the intensity resulting from thermal population 
of excited oscillator states. 

It is clear from Fig. [5] that no major features ap- 
pear or disappear upon heating to 300 K and 550 K. 
The corresponding changes are consistent with modest 
temperature-induced damping of the phonon oscillator 
modes, and perhaps with some softening of their en- 
ergy (ies). However, as shown by the FHWM ellipses in 
Figures 0] and El the resolution of the present measure- 
ment is too coarse for this issue to be carefully examined. 

To establish a quantitative detection limit for co- 
herently fluctuating local spin S e ff, we fit constant- E 
phonon scans through the longitudinal phonon near the 
(0,0,2) Bragg reflection to Gaussian profiles, Fig. [6l This 
yields phonon velocity v p h — 57.0(3) meV/rlu and the 
energy-integrated intensity of 280(50) cts-meV, which 
corresponds to the calculated cross-section of 110(5) 
mbarn (see Appendix lAl for details). Detection limit of 
100 cts-meV on the measured energy-integrated inten- 
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FIG. 6: Constant-energy scans through longitudinal phonon 
mode near (0,0,2) Bragg reflection at E = 8 and 12 meV at 
T = 6 K, which are part of the contour map shown in Fig. [5] 
(a). Curves are fits to Gaussian peaks. Small peak at I m 2.05 
in the E — 8 meV scan is a tail of the Bragg reflection. 



sity of magnetic mode, corresponding to twice the error 
bar obtained from phonon fits, gives detectable magnetic 
cross-section of > 39(2) mbarn. For |.F m (q)| 2 > 0.4, 
which holds for substantial part of our measured data, 
this leads to a detectable limit of S^j > 0.25. Thus, 
our results establish an upper limit on single-mode mag- 
netic scattering of less than 13% of the expected magnetic 
intensity for a singlet-triplet transition in the local-spin 
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FIG. 7: Partial density of states S P (E) obtained by numerical 
integration of the scattering intensities shown in Figure [4] 



picture. 

In order to quantify the temperature-induced changes 
more precisely, we have carried out numerical integration 
of our data at each energy with respect to the wave vector 
covered in each map, so as to obtain the partial scatter- 
ing density of states S P (E) shown in Fig. and the 
corresponding local dynamical imaginary susceptibility 
Xp(E) shown in Fig. [H In both cases we observe no in- 
dication of coherent peak corresponding to a well-defined 
localized-spin magnetic excitation. The error bars of our 
measurement constrain its possible intensity to < 20% of 
the corresponding phonon intensity. It should be men- 
tioned here that in related itinerant-electron pnictides 
where localized 3d magnetic moments are present, such 
as MnSb and CrSb, intense magnetic excitations on par 
with phonons were observed by neutron scattering.— £ 

Our data, however, do not rule out a possibility that 
a broad continuum of magnetic excitations of significant 
integral intensity exists, corresponding to delocalized cor- 
related magnetic states of itinerant electrons. In fact, a 
hint of such a continuum is contained in the tempera- 
ture evolution of Xp(E) m Fig- 13 which does seem to 
decrease slightly upon heating in the broad energy range 
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FIG. 8: Partial x'p(E) obtained by numerical integration of 
x"(Qi E) shown in Figure(5] Dashed lines are phenomenolog- 
ical fits to x"(E) of the damped harmonic oscillator. They 
primarily serve as guides for the eye, quantifying decrease of 
x"{Qi E) with temperature in the range 4 < E < 60 meV cov- 
ered in our measurement. Insets in (a)-(c) show T-dependence 
of the partial oscillator strength, 7 U = / 5m ^ ^ ' x"(E)dE 
calculated using data in main panels. 
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between 5 meV and 35 meV. Although this decrease is 
not clearly marked outside the statistical error of each 
given point of our measurement, there seems to be a 
consistent trend between different points, as well as for 
Xp{E) corresponding to maps at different Q shown in 
panels (a) - (c). It is also supported by phenomenolog- 
ical fits of our data to \" °f t ne damped harmonic os- 
cillator (DHO), shown by broken lines interpolating the 
data points in Fig. [H These fits are purely phenomeno- 
logical and are used simply to quantify the decrease of 
x"(Q,E) with temperature in the range 4 < E < 60 
meV covered in our measurement. Finally, it could be 
further quantified by the decrease of the partial (within 
our energy range) oscillator strength, I u , which is the 
quantity involved in the first moment sum rule for the 
dynamical spin susceptibility. 43 1^ calculated by numeri- 
cal integration of the corresponding data is shown in the 
insets in all three panels of the figure. The upper limit 
on the fluctuating itinerant spin could be estimated from 
the decrease in the integral intensity of magnetic scatter- 
ing of roughly < 100 counts- me V, which translates into 
S e ff ;$ 0.7, assuming the average magnetic form factor 
squared of ~ 0.2, Fig. [tJ] This is in agreement with 
the temperature- induced magnetic moment of 1.2/is ob- 
tained from static susceptibility measurement^ 

If indeed there is a weak continuum of magnetic scat- 
tering in the energy range between 5 meV and s» 35 
meV, its unambiguous identification is outside the lim- 
its of the presently available neutron technology. How- 
ever, it will become possible in the near future when new 
high-throughput polarized neutron inelastic spectrome- 
ters, such as HYSPEC at the Spallation Neutron Source 
at Oak ridge in the US, enter operation. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have conducted extensive search for 
magnetic scattering in the narrow-gap semiconductor 
FeSb2, which exhibits temperature- induced paramag- 
netism concomitant with the anisotropic increase of elec- 
trical conductivity. The temperature-induced paramag- 
netic susceptibility and the associated Shottky-like elec- 
tronic specific heat could be analyzed either in the frame- 
work of low (S=0) to high (S=l) spin state transition in 
the localized spin model^ or as a Pauli susceptibility 
of itinerant electrons in the itinerant narrow-gap picture 
similar to FeSi i 31 ' 37 A well-defined singlet-triplet mag- 
netic excitation between the crystal field split spin states 
is expected around ~ 50 meV in the localized spin pic- 
ture. Only a broad continuum of magnetic excitations 
corresponding to correlated magnetic states of excited 
itinerant electrons would be expected in the band pic- 
ture. In both cases magnetic intensity is expected to ex- 
hibit characteristic temperature dependence, decreasing 
upon heating, as excited states get thermally populated. 

In our data we find no evidence for a well-defined mag- 
netic excitation corresponding to transitions between the 



non-magnetic ground state and states of magnetic mul- 
tiplet in the localized spin picture. We find that peaks 
in the inelastic response of FeSb2 for energies up to 60 
meV and at temperatures s» 4.2 K, ss 300 K and « 550 
K are essentially consistent with the scattering by lattice 
phonon excitations. Our data establish a quantitative 
limit of S%ff < 0.25 on the fluctuating local spin. 

However, a broad magnetic scattering continuum in 
the 15 meV to 35 meV energy range is not ruled out by 
our data. In fact, a hint of such a continuum could be 
traced in the temperature dependence of the imaginary 
part of the local dynamical susceptibility shown in Fig. 
[51 which decreases upon heating. There should also be 
an accompanying weak quasielastic paramagnetic scat- 
tering similar to that observed in FeSi^ 3 - The putative 
magnetic intensity, however, is very weak and its unam- 
biguous identification by the means of polarized neutron 
scattering would have to be postponed until future de- 
velopments of the neutron scattering technology. 

Our findings make description of FeSb2 in terms of the 
localized Fe magnetic states unlikely and suggest that 
paramagnetic susceptibility of itinerant electrons is at the 
origin of the temperature-induced magnetism in FeSb2 . 
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Appendix A: Magnetic and phonon scattering 
cross-section 

In order to provide quantitative estimates, we com- 
pare magnetic scattering expected for a localized, non- 
dispersive triplet mode corresponding to a singlet-triplet 
transition at an iron site, to scattering cross-section for 
the longitudinal acoustic phonon near (0,0,2) Bragg re- 
flection in Fig. [4] (a). Magnetic scattering cross-section 
for a resonant non-dispersive magnetic mode is, 

A^| Mq) |a, (£ . A)i ( A1) 

where fej and kt are the incident and the scattered neu- 
tron wave vectors, A is the mode energy, N is the num- 
ber of unit cells in the sample (factor 2 accounts for two 
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FIG. 9: (a) - (c) Fe 4+ ionic magnetic from factor squared 
for maps shown in Figures [4] [5] (d) Calculated magnetic 
scattering cross-section for a non-dispersive, isotropic local 
triplet mode with S e ff = g/2y / S(S + 1) = 1. (e) Calculated 
scattering cross-section for a phonon in FeSb2 polarized along 
the wave vector transfer q at E — 8 meV, to be compared with 
panel (d). 



triplet transition, S e ff = g/V% ~ y/2, and magnetic 
scattering cross-section is twice that shown in Figure [9] 
(d), which was calculated for for S e ff — 1. 

For the long-wavelength acoustic phonon, the scatter- 
ing cross-section at T = is given by, 



dEdfl 



(A3) 



where e(q) is its energy, and the scattering length 
squared is, 



\bph (q) I 



2 h 2 q 2 cos 2 (3 , N|2 o nn q 2 



2M ceH e(q) 



l^(q) 



2.09 



M ceU e(q) 



(A4) 

In the latter, q is in A , M ce u in atomic mass units, 
e(q) in meV, and /3 is an angle between the phonon po- 
larization vector and the wave vector transfer, q. F(q) 
is the unit cell structure factor, 



Fe ions per unit cell), and magnetic scattering length 
squared is given by, 

|6 m (q)| 2 = r 2 i |F m (q)| 2 ^ e 2 // . (A2) 

Here F m (q) is magnetic form factor for the corresponding 
magnetic ion (we use |F m (q)| 2 for Fe 4+ for our estimates, 
Fig. [9](a)-(c)), r m = — 5.39- 1CU 13 cm, and we have intro- 
duced the effective spin through Seff = g/2y/S(S+l), 
where g is the spectroscopic Lande factor. For a singlet- 



F(q) = ^2b ll e- i ^\ (A5) 

where fj, indexes atom at a position in the unit cell, 6 M 
is its nuclear scattering length. Scattering cross-section 
for the longitudinally (|j q) polarized acoustic phonon at 
E = 8 meV, calculated from Eqs. (|A3|) - (|A5I) . is shown 
in Figure [9] (e). It is already clear from comparison of 
Fig. IH1 (d) and (e) that in the range where our measure- 
ments were performed the expected intensity of the local 
magnetic mode in the local-spin picture is similar to, or 
higher than that of the longitudinal phonon in Fig. |4j 
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